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;h ! Abstract 
Ph. 

<^ I We study the C and P even WW'j and WW Z trilinear gauge boson vertices (TGV's), in 

(3 ■ the context of the MSSM as functions of the soft SUSY breaking parameters Aq, mo, M1/2 

^ ■ and the momentum q carried by 7, Z, assuming the external W's are on their mass shell. 

^ . We follow a complete renormalization group analysis taking into account all constraints 

^ ! imposed by the radiative breaking of the electroweak symmetry. It is found that for 
energies ^/s = < 200 GeV squark and slepton contributions to the aforementioned 



' couplings are two orders of magnitude smaller than those of the Standard Model (SM). 



In the same energy range the bulk of the supersymmetric Higgs corrections to the TGV's 
irS ■ is due to the lightest neutral Higgs, /iq, whose contribution is like that of a Standard 
^ ! Model Higgs of the same mass. The rest have negligible effect due to their heaviness. 
\ The contributions of the Neutralinos and Charginos are sensitive to the input value for 
I I the soft gaugino mass M1/2 being more pronounced for values M1/2 < 100 GeV . In this 
case and in the unphysical region, < ^/s < 2Mw their contributions are substantially 
^ . enhanced resulting to large corrections to the static quantities of the W boson. However 
such an enhancement is not observed in the physical region and their corrections to the 
' TGV's are rather small. In general for 2M\y < y/s < 200 GeV the MSSM predictions 
^ I differ from those of the SM but they are of the same order of magnitude. Deviations from 
the SM predictions to be detectable require sensitivities reaching the per mille level and 
hence unlikely to be observed at LEP200. For higher energies SM and MSSM predictions 
exhibit a fast fall off behaviour, in accord with unitarity requirements, getting smaller 
by almost an order of magnitude already at energies ^/s ~ .5 TeV. At these energies the 
task of observing deviations from the SM predictions, which are due to supersymmetry, 
becomes even harder requiring higher accuracies. 
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1. Introduction 



The Standard Model (SM) as been remarkably succeful in describing particles interac- 
tions at energies around the ~ lOOGeV. Precise measurernenis at LEP provided accurate 
tests of the standard theory of electroweak interactions Ef but we are still lacking a 
direct experimental confirmation of the non-abelian structure of the standard theory. 
The WlVj, WWZ and ZZj couplings are uniqely determined within the context of the 
SM and such couplings will be probed in the near future with high accuracy. The study 
of the trilinear gauge bosons vertices (TGV's) in the e"'"e~ — > W~^W~ process is the pri- 
mary motivation for the upgrading of LEP200 S and the potential for measuring these 
has been discussed in detail B 0. At an energy of about 190 GeV and with integrated 
luminosity of 500 an accuracy of .1 for the determination of these couplings can be 
obtained. So far there are no stringent experimental bounds on these couplings and 
the efforts of the various experimental groups towards this direction are still continuing. 
In the near or remote future with the proposed or already under construction high energy 
colliders (LHC, NLC, CLIC, JLC) further improvements on the TGV's bounds will be 
obtained reaching accuracies Such a precise measurements are of 

vital importance not only for the SM itself but also for probing new physics which opens 
at scales larger than the Fermi scale. 

The gauge boson vertex has been the subject of an intense theoretical study the last 
years. In particular the WWV vertex (1/ = 7 or Z) has been analysed in detail within 
the framework of the standard theory, as well as in extension of it, and its phenomenology 
has been discussed. The lagrangian density describing the WWV interaction is given by 

Cwwv ^ _,g^^^^^^rl^^r,v^ _ + KyW^W^F^^ + ^wl^Wj^V^^ + ...], (1) 

f e for y = 7 

where the ellipsis stand for P or C odd terms and higher dimensional operators. In 
Eq. (1) the scalar components for all gauge bosons involved have been omitted, that is 
d-W = d-V = 0, since essentially they couple to massless fermions []. At the tree level 
Kv and Xy have the values /ty = 1, Ay = 0. However radiative corrections modify these, 
the order of magnitude of these corrections being 0{^) ~ 10^^. Sensitivity limits of this 
order of magnitude will not be reached at LEP200 but can be achieved in future colliders 
where the TGV's can be studied in detail and yield valuable information not only for 
the self consistency of the SM but also for probing underlying new physics. Any new 
dynamics whose onset lies in the TeV range modifies Ky, Ay and deviations from the 
SM predictions are expected. 

Supersymmetry (SUSY), is an extension of the SM which is theoretical motivated 
but without any experimental confirmation. The only experimental hint for its exis- 

^For on shell VF's we have (□ + m'^)W^ = and certainly d-W — 0, while for the photon d-A 
vanishes on account of current conservation. 
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tence derives from the fact that the gauge coupHngs unify at energies ~ 10^^ GeV if 
we adopt a supersymmetric extension of the SM in which SUSY is broken at energies 
MsusY^O{l)TeV . Supersymmetric particles with such large masses can be produced 
in the laboratory provided we have very high energies and luminocities. However their 
existence affects Ky and Ay, even at energies lower than the SUSY production threshold 
making them deviate from the SM predictions. Therefore the study of these quantities 
may furnish as a good laboratory to look for signal of supersymmetry at energies below 
the SUSY production threshold. In any case such studies serve as a complementary test 
along with other efforts towards searching for signals of new physics and supersymmetry 
is among the prominent candidates. 

In the SM Ky and Ay as functions of the momentun carried by the V boson 
{V = 7,2'), for on shell W^s, have been studied in detail but a similar analysis has 
not been carried out within the context of the MSSM. Only the quantities = 0), 

A^(g^ = 0) have been considered which are actually related to the static quantities 
magnetic dipole (fiw) and electric quadrupole (Qw) moments of the W boson. To 
be of relevance for future collider experiments the form factors n-y^z, ^■y,z should be 
evaluated in the region > 4M^. The behavior of k^^z, \,z in this physical region 
may be different from that at = especially when the energy gets closer to Msusy 
and supersymmetric particles may yield sizable effects, due to the fact that we are 
approaching their thresholds. In those cases an enhancement of their corresponding 
contributions is expected, unlike SM contributions which in this high energy regime are 
suppressed. We should also point out that some of the supersjTumetric particles may have 
relatively small masses, for a certain range of the parameters and in those circumstances 
their contributions to TGV's are not necessarily small. In order to know the magnitude 
of these effects a detailed computation of the trilinear gauge boson couplings should be 
carried out. 

In this work we undertake this problem and study the C and P even WW'y, WW Z 
vertices in the context of the MSSM when the external W bosons are on their mass 
shell. Such studies are important in view of forthcoming experiments at LEP200 and 
other future collider experiments which will probe the structure of the gauge boson 
couplings and test with high accuracy the predictions of the SM. If deviations from the 
SM predictions are observed these experiments will signal the presence of new underlying 
dynamics which opens at scales larger than the Fermi scale. 

The magnitude of the aforementioned couplings and their dependence on the ar- 
bitrary parameters of the MSSM requires a systematic study in which all limitations 
imposed by a renormalization group analysis of all running parameters involved and es- 
pecially those arising from the radiative breaking of the electroweak (EW) symmetry are 
duly taken into account. In this paper we deal with this issue and calculate the trilinear 
vector boson couplings as functions of the momentum carried by the 7, Z and the 
arbitrary parameters of minimal supersymmetry. 



2 



This paper is organized as follows: 
In section 2 we give a brief outline of the MSSM. In section 3 we carry on to discuss the 
SM predictions for the TGV's paying special attention to the contributions of fermions 
and those of the gauge bosons discussing the issue of gauge independence. In section 4 
the MSSM predictions for these couplings is discussed. Section 5 deals with the absortive 
(Imaginary) parts of these vertices and in section 6 we present the numerical analysis 
and discuss our conclusions. 



2. The MSSM 

The MSSM is described by a Lagrangian 

1^ = ^SUSY + ^soft (2) 

where its supersymmetric part Csusy is derived from a superpotential W bearing the 
form 

W = ihuQ'MU' + hDQ'MD' + hEL'ME'' + fiHiHihj , = +1 (3) 

In Eq. (3) carrets denote supermultiplets. Minimality is enforced by assuming the SU (3) x 
SU{2) X f/(l) as the gauge group and the least number of chiral multiplets necessary 
to accomodate matter fermions and drive EW symmetry breaking. The superpotential 
above conserves i?-parity. The Hi, H2 Higgses give mass to up and down fermions re- 
spectively after EW breaking takes place. The part responsible for the soft breaking of 
supersymmetry is given by 

-Csoft = Y.^im'' + (huAuQH2W + hDADQHiD' + hEALLHiE' + h.c.) 

i 

+ {fiBHiH2 + h.C.) + \Y,Ma\a\a- (4) 

where rrii are the soft scalar masses, Au,d,l are trilinear soft couplings, B is the Higgs 
mixing parameter and , a = 1, 2, 3 are the soft gauginos masses for the f/(l), SU{2) 
and S'f/(3) gauge fermions respectively. Throughtout this paper we assume universal 
boundary conditions at the unification scale Mqut — 10^^ GeV 

ml = ml , Au^D^L ^Aq , M„ = M1/2 (5) 

This choice is suggested by grand unification and also by absence of flavor changing 
neutral currents which puts stringent constraints on the difference of masses squared of 
the same charge squarks. However this is in no way mandatory. At least it parametrizes 
our ignorance concerning the origin of the soft SUSY breaking parameters in the most 
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economical and plausible way. We don't expect that altering the boundary conditions 
given in Eqs.(5) above will drastically affect the estimates for the TGV's. 

The values of all running parameters in the vicinity of the electroweak scale are 
then given by solving their Renormalization Group Equations (RGB's) having as initial 
conditions Eqs. (5) 

The breaking of the EW symmetry is known to proceed via radiative corrections 



driven by the large top Yukawa coupling 1121. The equations minimizing the scalar 
potential of the theory are 

— = tan' (3 _ 1 ' '''''^f^ = ~m\+ml 

where the angle (3 sets the relative strength of the v.e.v. of the Hi and H2 Higgs fields 
involved 

In Eqs. (6)-(7) all quantities are meant at M^, the experimental value of the Z boson 
mass {Mz = 92.18 GeV), and , are defined as 

- 2 2 I dA.V 2 2,2 /o^ 

"^1,2 = "^1,2 + 5 rn^,2 = ^H,2+ ^^ ■ (S) 

OV12 

In Eq.(8) mHi2 s-^'^ the soft Higgses masses and mi 2 differ from mi^2 by the contributions 
of the one loop corrections AV to the scalar potential of the theory. Including the one 
loop corrections within the minimizing Eqs. (6), as prescribed by Eq.(8), is a necesssary 
ingredient for the numerical stability of our physical results. If not included the physical 
quantities would strongly depend on the choice of the scale at which physical quantities 
are evaluated leading to results that are ambiguous and untrustworthy 'c^. 

The arbitrary parameters of the model are the soft parameter mo, A, M1/2, B the 
mixing parameter fi as well as the values of the top Yukawa coupling ht 0. This number 
is reduced to five if use is made of the first of the minimizing equations (6). Then a 
convenient set of independent parameters, which is adopted by many authors, is to take 
the set 

mo , Ao , Ml/2 , tanp{Mz) , mt{Mz) (9) 

where mt{Mz) is the value of the "runing" top quark mass at the scale Mz- This 
facilitates the numerical analysis a great deal since the RGE's of all soft masses and 
parameters involved, with the exception of B and /i, do not depend on B, /i (nearly 
decouple). Given the inputs (9) the RGE's can be solved and predictions for the mass 
spectrum can be given. At this point we should remark that all subtleties associated 
with this approach, like for instance the presence of low energy thresholds and other 



^The top Yukawa coupling although localized in a narrow range of values in view of the recent exper- 
imental evidence for the top quark 113 is being considered as an input parameter since the occurrence 
of symmetry breaking is very sensitive to its value. 
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uncertainties due to higher loop effects, in no way affect the one loop corrections to the 
TGV's. Complete expressions for the RGE's of all parameters involved can be traced in 
the literature and will not be repeated here 0' Q. 

After this brief outline of the MSSM we embark to discuss the TGV's defined in the 
previous section. 



3. The SM contribution to Akv{Q^), AQv(Q^' 

Altho ugh the SM contributions to the TGV's have already been calculated in the liter- 
ature tcJ, for reasons of completeness we shall briefly discuss them in this section too 
paying special attention to the contributions of fermions and gauge bosons. 

In momentum space the most general WWV vertex (V = 7 or Z) with the two W's 
on shell and keeping only the transvers degrees of freedom for the 7 or Z can be writen 
as0 

+2 Akv {ga^^Qp - gp^^Qa) + 4 A^(Q,g^ - ^g^p)} + ... (10) 

where Aky = ky + ~ ^ , ^Qv = ~2 Xy- The kinematics of the vertex is shown in 
Fig. 1 . 

The ordinary matter fermion contributions both to = and ^ have been 
studied elsewhere IHI HH . However in those works there is an important sign error which 
affects substantially the results given in those references This has been also pointed 
out independently in ref. . The consequences of this for the static quantities of the 
W boson fiyr, Qy/ has been discussed in detail in ref. [|17|, 0. In the massless fermion 



limit, which is actually the case for the first two families, this leads to nonvanishing 
contributions for both Ak^lQ"^ = 0) and AQ^IQ"^ = 0), contrary to what had been 
previously claimed. These are proportional to Tr{QT^) 7^ 0, unlike the anomaly terms 
which are proportional to TrQ and hence vanishing. The details of this calculation 
which points out this important sign error can be traced in the literature In the 
SM these contributions to Ak^{Q'^ = 0) are large and negative partially cancelling the 
contributions of the gauge bosons and the standard model Higgs which are positive. 

In units of g'^/lQir'^ the fermion contributions of the triangle graph shown in Fig. 2a 
are as follows, 

Akv = -cvTiCg j^dt j^da{gi[t^ + t^{-l + Rf: -Rf)+t^{Rf-Rfi) 



AQv - -cyTiC,gifdtfda '''^'-f-^^^ (12) 

Jo Jo Lj 
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where c^ = l,cz = R with R = M^M^. In Eqs. (11) and (12) 



Lj^t'{l-^a{l-a))+tiRf-Rf>-l) + Rf,+te , RfJ' ^ 



and Cg is the color factor (1 for leptons, 3 for quarks). The couphngs gl appearing in 
Eq.(ll) are as follows 



V = 1 ■ 9l = 9r = QL 

V = Z : gi = Qi^,gf^ = Ql^ 

where Qem are the electromagnetic charges, and r are the weak charges for left /right 
handed fermions, which are defined by the relation = — Qf^sin^ 6w is the 
weak isospin of the fermion /, ie. = —1/2 for the left handed charged leptons etc. At 

2 

zero momentum transfer Ak^, AQ^ are nonvanishing as said earlier yielding -^^0{1) 
contributions to both A/c^, AQ^. Actually it turns out that when = this is the larger 
contribution to AQ^ of all sectors, while Ak^{Q'^ = 0) is negative partially cancelling 

2 

large (~ -j^O^l)) contribution from gauge boson and Higgs particles as said previously. 
We should point out that Aky, AQy as given in Eqs. (11) and (12) refer to both the real 
and imaginary parts of the vertices; as we pass the internal particle threshold imaginary 
parts develop, too due to the ie appearing in the denominator. One can easily check 
that as — ^ cxD fermion contributions to both Aky, AQv, tend to zero as demanded 
by unitarity. 

Regarding the contributions of the gauge bosons to Aky, AQy the calculations were 



carried out in the 't Hooft - Feynman gauge IIj , and the results are known to be gauge 
dependent. The details of this calculation can be traced in the literature [see Eq.(8)- 
(23) of ref. 0]]. In order to render the trilinear gauge boson vertices gauge independent 



we should add to them additional contributions from box graphs by applying special 
field theory techn iqu es, such as the pinch technique 0, or work in manifestly gauge 
invariant gauges 1^. As a result of this gauge dependence the quantity Aky turns 
out to have bad high energy behaviour, growing logarithmically as the energy increases, 
violating unitarity constraints. Besides being gauge dependent Aky is also singular at 
the infrared [IR). Actually this IR singularity occuring in one of the graphs is the only 
one surviving among several other which cancel against each other. As said previously 
for the restoration of gauge independence additional contributions from some box graphs, 
the pinch contributions, should be appended to the vertex parts. These also cancel the 
IR divergence mentioned earlier. In units of gf^/lGyr^ these pinch parts are given by 'S, 
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At 



1 ^2 
t2 + _ ^) _ 



2t 



-t'^ail — a) + it 



+ sin^ 6*^^ / (it— 
Jo 1 



-t(l -t) + ze 



(13) 



1 

- R){cosHw I dt 



da- 



2t 



+ R{1 -t)- ^t^a{l -a) + ie 



+ sin^ Ow j dt— 
Jo 1 



1 



'-t{l-t)+ie 



} + IR 



(14) 



In Eqs. (13)-(14), R = Ml/M^ while IR stands for the infrared singularities which 



cancel against the corresponding singularities of the vertices given in ref. [n]. Once the 
pinch contribution Eq.(13)-(14) are taken into account the gauge boson contributions 
become gauge independent approaching zero values as increases as demanded by 
unitarity and are also free of infrared singularities. 



4. The MSSM contribution to A/cv(Q2), AQy(Q2) 

At = the MSSM contributions to A/c^, /S.Q^ have been studied elsewhere and their 
dependences on the soft breaking parameters A, niQ, M1/2, tan/5 and top quark mass mt 
have been investigated The q, I (squarks, sleptons), Z, C (neutralinos, charginos) 
as well as the supersymmetric Higgs contributions to Aky, AQv are deduced from the 
triangle graphs shown in Figures 2 and 3. We show only graphs that yield nonvanishing 
contributions to at least one of the A/cy, AQy. In the following we shall consider the 
contributions of each sector separately. 

Squarks— Sleptons {q, I) 

We first consider the contributions of the sfermion sector of the theory which can be read 
from the diagram (b) of Figure 2. Unlike matter fermions this graph involves mixing 
matrices due to the fact that left fi and right /£ handed sfermions mix when electroweak 
symmetry breaks down. Such mixings are substantial in the stops, due to the heaviness 
of the top quark, resulting to large mass splitting of the corresponding mass eigenstates 
ti^2- For the sfermions we find after a straightforward calculation that: 
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Ah 



dt 



1 e{i-t){2t-i 

da 



Ri 



■I i 







dt 



4g^ 



1 , 4t3(l-t)a(l 
da 



7'2 



a 



^2 



(15) 
(16) 



Lj = r + {Rf -Rp-l)t + Rp~-^ra{l-a)+ie , Rf. ?, = {mf. ?, /Mw] 



w 

The prefactors appearing in the integrals above cy,T3 , g]^ and Cg, are exactly those of 
fermions, see Eqs. (11)-(12), since fermions and their superpartners carry same quan- 
tum numbers. /i,2 and /{ 2 denote the mass eigenstates while K-'^'-'^' diagonalize the 

corresponding mass matrices, i.e. K-^'-^ A1^^,K-^'''^ = diagonal. In the stop sector for 
instance, where such mixings are large, the corresponding mass matrix is given by 



Mi 



+M|(cos2/3)(i - lsin^9w) 

V 



mt{A + ncot[3) 



\ 



(17) 



+Ml{cos2l3){\sm'^ew) j 



and the diagonalizing matrix is defined as K*A^?K* = diagonal{ml^^ml^). In the 

absence of SUSY breaking effects mj j, = rajj' and K-^'-^' become the unit matrices. In 
that limit AQy given above cancels against the corresponding fermionic contribution as 
it should. 



Neutralinos— Charginos (Z, C) 

The neutralino and chargino sector is perhaps the most awkward sector to deal with 
owing to mixings originating from the electroweak symmetry breaking effects. Their 
contributions are read from the graphs shown in Figure 2c, d. In the following we shall 
denote by Ci the two chargino states (Dirac fermions) and by the four neutrali- 
nos states (Majorana fermions). Recall that they are eigenstates of the following mass 
matrices 113, 

■Mc=f ~'"'] (18) 

^ \ -gvi 
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(19) 



/Ml g'vi/V2 -g'v2/V2 \ 

M2 -gvi/V2 gv2/V2 

g'vi/V2 -gvi/V2 

\ -fifS/Vs gv2/V^ -/jl y 

where fi/v/2, V2/V^ arc the v.e.v.'s of the neutral components of the Higgs field i^i and 
H2 respectively. If the U, V matrices diagonalize A^g' i-^- UA^qV^ = diagonal, and 
O (real orthogonal) diagonalizes the real symmetric neutralino mass matrix of Eq.(19), 
O^A4f^O = diagonal, then the electromagnetic and weak currents are given by 



i 



(20) 

(21) 
(22) 



where 



riR 



B. 



L 

a/3 



— --^03aU*2 - 02aU*i , C^j — +-j=04aV^ - 02aV*i 

= [cos^ 9wSij - ^{Vi2V*2ShL + U*2Uj2ShR)] , h = L,R 



(23) 
(24) 
(25) 



Pb..l arc the right/left handed projection operator (1 ±75)72. The contributions of this 
sector to Ak^, /^Q^, as calculated from the graph shown in Figure 2c, is as follows: 



-T. I' dt [' da{Fai[t^ + {Ra - - l)t^ + {2Ri - R^)^ 
10 Jo 



^ t\3t-2)a{l-a)] + Go 



f^ima, 2 - - 1 



'-{At'-2t)}- 



-S^F^, I d.t I da- ^ ^ ^ ' 



JO 



4 



(26) 
(27) 



r + {Ri-Rc,- l)t + Ra- jj^ra{l -a)+i€ 



C. 



R |2 



C. 



L |2 



iC^fi^ + h.c.) 



(28) 
(29) 



where the index a = 1, 2, 3, 4 is the neutralino index and i = 1, 2 is the chargino index. 
Note that we have not commited ourselves to a particular sign convention for the masses 
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rrii, rua appearing to the sum in the equation above for the Ak^. Chiral rotations that 
makes these masses positive it also affects the rotation matrices and should be taken into 
account. 

For the couplings Akz, AQz we get, 
Graph of Figure 2c : 

Akz = -R f dtj^ da{S('j^[f{l-t){-t + Ria + Rj{l-a)-Ra) 

-'^iT,% + T,'}Jl2t'a + e-t]}-±- (30) 

nl rl ^3 



where 



AQ. = -8fl E S^, /' * /' ,„ t^(l-'Hl-«) (31) 



402 

Lj^ = atRi + t{l - a)Rj + - t) - t{l - t) - jj^t^a{l - a) + 



le 



2 

m 



Rij^a = '2 , i,j — chargino indices , a — neutralino index 



Graph of Figure 2d : 

Same as in previous graph with a} replaced by {p, cr, i} and Sl^j^^\T^p^\ Lfj^ re- 
placed by the following expressions: 



S''f^ ^ -{Ci:C^,Bi^^^ + {L^R)) 

T'T ^ -iC5*C^iBi(^^ + {L^R)) 

Ljf ) = atRp + t{l - a)Ra + Ri{l - t) - t{l - t) - -^t^a{l -oc)+ie 

a, p = neutralino indices , i = chargino index 



Higgses {Ho,ho,A,H±) 

There arc five physical Higgs bosons which survive electroweak symmetry breaking. Two 
of these, Hq and ho, are neutral and CP even, while a third A, is neutral and CP 
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odd. The remaining Higgses bosons, H±, are charged. At the tree level the lightest of 
these, namely Hq, is lighter than the Z gauge boson itself. However it is well known 
that radiative corrections which are due to the heavy top are quite large and should 
be taken into account. These modify its tree level mass by large amounts 5m\^ ~ 
gi'(mj/M^) ln(m|^m|^/mj) which may push its mass up to values exceeding Mz ; in 
some cases up to ~ 130 GeV. ho turns out to yield the largest contributions of all 
Higgses to the TGV's since the remaining Higgses have large masses of the order of the 
SUSY breaking scale. At the tree level the masses of all Higgs bosons involved are given 
by the following expressions : 

m\ = , {ml^B^^) (32) 



„2 



sin 2/3 

^{{ml + Mlf ± ^J{m\ + M^f - AMlm\cos^{2P)} (33) 
rriH^ - m\ + M^ (34) 

The Higgs contributions can be expressed in terms of their masses and an angle 9, which 
relates the states 5*1 = cos j3 [Real H^) + sin j3 {Real H2) , 5*2 = — sin (3 [Real Hf) + 
COS/3 {Real H^) to the mass eigenstates ho,Ho. The state Si is the SM Higgs boson 
which however is not a mass eigenstate since it mixes with 82- When sin' 9 — 1 such a 
mixing does not occur and ho becomes the standard model Higgs bosons 5*1. 

The contributions of the Higgs bosons to Ak^, AQ^, follow from the graphs shown 
in Figure 3 and are as follows. 
Graphs of Figures 3a and 3b : 



A: Ak^ = D2{Ra, R+) , AQ^ = Q{Ra,R+) (35) 

ho : Ak^ = sm^9Di{Rh„) + cos" 9 D2{Rho, R+) (36) 

AQ^ = smHQ{Rh„l) + cosH Q{Rh„R+) (37) 

Hq : As in ho with Rh^ — > Rh^ and sin' 9 ^ cos' 9 (38) 

T, -f /M h TJ A n -2. Ml + M| Sin' 2/3 -M,'„ 
Ra = {ma Mw) , a = ho, Ho, A, H± , sm 9 ^ —2 . 

The functions D12, Q appearing above are defined in Appendix A. 

For the Akz, AQz form factors there are more graph contributing. These are shown 
in Figures 3c-3f. From the graphs shown in Figure 3 we pick the following contributions. 
Graph of Figure 3a: 

Akz = i{(cos' 9) C dt C da [(4 - 2i?)t^ + {R- 2){Rh, + 2)t^ 
4 JO Jo 

+{2Rh, - RRho + 8 - 2R)t^]-^ + (sin' 9) x {Ho ho)} (39) 
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AQz = (2 - i^){(cos^ 9) [' dt [' d J'^^ ^j?'^^ "'^ 

Jo Jo Lir^ 



'Ho 

,2 



+isin'e) X (Ho^ho)} (40) 



Graph of Figure 3b : 

Akz = {^^){D2(Ra, R+) + sin' e D2(Rho, R+) + cos' D2(Rho,R+)} (41) 
AQz = {'^^){Q{Ra, R+) + sin' 6 Q{Rho, R+) + cos' 6 Q{Rho, R+)} (42) 

Graphs of Figures 3c and 3d : 

■ 2n^ ^ ^'(1 - ^)(1 + R+- RAOi - RHo(i - a) - 2t) 



7"2 



Akz = -{(sin^^) [ dt f da 
2 70 ^0 

+(cos'^) X (i/o^/io)} (43) 

AQz = 2i?{(sin2 6) C dt C da "^^^ ~ ~ + (cos' ^) x (i/o ^ /io)} (44) 

^0 Jo Li, 



40^ 

L^^ = -t{l-t) + RAat + RHo{'^-a)t + R+{l-t)-^^t^a{l-a)+ie{Ab) 



and finally, 

Graphs of Figures 3e and 3f : 



R 

A/cz = —{(cos^e) dt da[-6at^ + (t^-t^)(2(t-l) + (R-RHo)a + RHo) 
2 JO Jo 

+2{R - l)at']^ + (sin' 9) x (i/o ^ /^o)} (46) 



AQz = 2i?{(cos' 9) r dt [' da^ — — ^} + (sin' 9) x {Hq ^ ho)} (47) 

° JO L|^^ 



,2 



40 

= {1 - t)"^ + Rat + RHot{l - a) - -^t^a{l - a) + ie (48) 



In the most of the parameter space the Higgses A, H± and Hq turn out to be rather 
heavy having masses of the order of the SUSY breaking scale; therefore all graphs in 
which at least one of these participates are small. At the same time sin^ 9 has a value 
very close to unity. Thus the dominant Higgs contribution arises solely from the graphs 
of Figures 3a,e and 3d in which a ho is exchanged. This is exactly what one gets in the 
SM with ho playing the role of the SM Higgs boson. 
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5. The Absortive Parts of the TGV's 



The contributions of the TGV's presented so far have also imaginary (absortive) parts 
which show up as soon as one passes the thresholds associated with the particles ex- 
changed in the one loops graphs. Since the majority of the one loop expressions en- 
countered have a triangle structure some of these thresholds can be anomalous and this 
depends on the masses of the particles circulating in the loop. 

The absortive parts can be readily calculated using the ie prescription. Actually the 
denominators of the Feynman integrals involved carry a small positive imaginary part 
and it is a matter of a proper algebraic manipulation to pick up the relevant imaginary 
parts of all integral expressions presented in the previous sections. All graphs yielding 
nonvanishing contributions to Ak^^z, AQ^^z have a triangular structure with the excep- 
tion of some of the gauge boson graphs which involve the quartic gauge boson coupling. 
For these we have AQ^ z = with Ak^^z given by 



Ak 



7,^ 



ysm 9w 



+ 



t dt ln(l - 
Jo 



3 ^2t^ - + R)t^ + mt 

2RJo ^ t2 + jif^i -t)+ie 

2 



4Q- 



t a(l — a) + ie) 



(49) 



where [R = M'^/M^). The first of the integrals in the expression above does not have 
any discontinuity since the denominator never vanishes f\ The second integral developes 
an absortive part when the argument of the logarithm becomes negative. Therefore it 
has an Imaginary part given by 



ImAk^,z = —Q{Q'-M^] 



\ 



'1 



(50) 



For the pinch contributions we have an absortive part arising from the single dt in- 
tegrations appearing in the Eqs. (13), (14) and an additional contribution which stems 
from the double dt, da integrations of these equations. The later yield absortive parts, 
denoted by A^^z, which have a structure akin to those of the triangle graphs. The former 
yield absortive parts which are easily calculated leading to the following results: 



Im Ak^ 



ImAkz 



. 2, e(Q^-M^) 

vr sm t^i^ — , + 



TT sin^6' 



w 



R 



9(g^ 

4g2^ 



m2 



W) 



(51) 
(52) 



^ Actually this integral arises from a diagram that has the structure of a two point Greens function. 
It involves the quartic gauge coupling where one of its legs is the 7 or Z and the other is one of the on 
shell external W's. Since the W^s are on their mass shell it does not depend on Q^. 
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The absortive parts of all triangle graphs as well as the contributions lS.^^z in the 
Equations (51) and (52) above can be inferred by calculating the imaginary part of the 
integral 



I = C dtl 



1 



da 



+ Pia + P] 



(53) 



JO JO 



p ^ aa^ j^t'^a{l — a) + ie 



w 



where Po,i,2 as well as p, a are functions of the variable t alone. The absortive part of 
this integral is not difficult to calculate, and the details are presented in the Appendix 
B. In its final form it can be expressed as an single integral over the variable t (see Eq. 
(58), Appendix B) which can be integrated numerically using special numerical routines. 
Having expressed the Imaginary parts of all contributions involved, as integrals in t of 
known functions of t and the energy variable Q^, we are ready to proceed to numerical 
computations. 

The strategy we follow for the evaluation of the absortive parts will also apply to 
the real (dispersive) parts of the form factors under consideration. In fact wherenever 
a double Jq dt Jq da integrations are encountered we first perform the da integrations 
explicitly and subsequently evaluate the dt integrations numerically for various values 
of and the MSSM parameters. The details of the numerical analysis is the subject of 
the following section. 



6. Numerical Analysis — Conclusions 

As discussed in the previous section both dispersive and absortive parts of the trilinear 
WW J, WWZ vertices can be cast as single integrals of known functions of t and Q^, 
which also depend on the physical masses of all particles involved. These integrations 
we have numerically carried out using special routines of the FORTRAN Library IMSL 
available to us. The advantage of using this facility is that it leads to reliable results 
even in cases where the integrands exhibit fast growth at some points or have a rapid 
oscillatory behaviour. The inputs in these calculations are the value of the energy variable 
Q"^ and the arbitrary parameters of the MSSM discussed in the section 2. 

In our numerical analysis we have taken all Yukawa couplings, but that of the top 
quark vanishing, which is a very good approximation especially for the reason that the 
trilinear gauge boson vertices under consideration are already of one loop order. This 
approximation however holds provided the value of the parameter iaxi l3{Mz) which sets 
the relative strentgth of the v.e.v's of the two Higgses involved are not large < 10. For 
larger values the bottom Yukawa coupling should be also considered in the RGE's of all 
running parameters involved. However this approximation little affects our numerical 
results for the form factors under consideration. 

With the experimental inputs Mz = 91.18 GeV, sin^ Ow = -239, a^miMz) = 1/129 
and as{Mz) — .117 and with given values for the arbitrary parameters tan (3{Mz), 
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■mt{Mz), Ao, mo, M1/2 we run our numerical routines in order to know the mass spec- 
trum and the relevant mixing parameters necessary for the evaluation of the form factors 
given in the previous sections. Throughout the analysis we have taken tan/5(M^) < 10, 
but its value cannot be taken arbitrarilly small. Actually given mt{Mz) the parameter 
ia.n (3{Mz) is forced to a minimum value otherise Landau poles are encountered mak- 
ing the top Yukawa coupling /ij/(47r)^ getting values > 1 outside the validity of the 
perturbative regime. 

For the running top quark mass mt{Mz) we took values in the whole range from 
130 GeV to 190 GeV, although small values of rrit are already ruled out experimentally 
in view of the recent CDF and DO results which both quote a mass for the top quark 



larger than about 170 GeV 0. The physical top quark masses emerging out are slightly 
larger by about 3% Q . 

As for the soft SUSY breaking parameters Aq, mo, M1/2 we scan the three dimensional 
parameter space from ~ 100 GeV to 1 TeV. This parameter space can be divided into 
three main regions: 



i) Aq ~ mo — Ml/2 (all SUSY breaking terms comparable) 

ii) Aq ~ mo <^ Ml/2 (the gaugino mass is the dominant source of SUSY breaking ) 

iii) Ml/2 <^ ^0 — ""^0 (^o^'^o dominate over M1/2) 

Case ii) covers the no-scale models for which in most of the cases the preferable val- 
ues are Aq = uiq = while case iii) the light gluino case. 

Regarding the values scanned for the energy variable we moved both in the timelike 
and spacelike region for values ranging from | | = to \J\ Q"^ \ = lO^Myy. For the 
timelike case, which is of relevance for future collider experiments, this corresponds to 
values of a/s ranging from GeV to about 600 Mw- For comparison we quote that -y/i 
at LEP200 will be 190 GeV that is it just exceeds the two W^s production threshold 
energy 2Mw- Both in the spacelike and timelike energy region as soon as ^/s exceeds 
~ few TeV the contributions of each sector separately becomes negligible approaching 
zero as the energy increases in accord with unitarity requirements. 

Sample results are presented in Tables I and II for values of (Ao, mo, M1/2) equal to 
(300, 300, 300), (0, 0, 300) and (300, 300, 80) GeV representative of the cases i),ii) and iii) 
respectively discussed above. The inputs for the remaining parameters are tanjd^Mz) = 
2, mt{Mz) = nOGeV. The value of ^/s in these tables are respectively 190 and SOOGeV^, 
corresponding to the center of mass energy of LEP200 and NLC In the same tables 
for comparison we give the SM predictions for Standard Model Higgs masses 50, 100 
and 300 GeV. With the inputs given above the typical SUSY breaking scale lies in 
somewhere between 2Mw and .5 TeV. Although many sparticle thresholds exist in this 
region, as for instance the lightest of the sleptons and squarks as well as the lightest 
of the neutralinos and charginos, especially when M1/2 is light, these thresholds do not 

Our results for Aky, AQy are not sensitive to the choice of the angle (3. 
^ The physical top quark mass Mt in the DR scheme is given by Mt = mt{Mt) / [1 + 5a3(Afi)/37r). 
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result in any enhancement of the form factors Ak^^z, AQ^^z- Increasing the value of 
the dominant SUSY breaking scale the supersymmetric contributions to these quantities 
become less important approaching zero values. Of all sectors the sfermions yield the 
smaller contributions in the entire parameter space even in cases where due to large 
electroweak mixings some of the squarks, namely one of the stops, are relatively light. 
The supersymmetric Higgses yield contributions comparable to those of the SM, provided 
the latter involves a light Higgs with mass around the 100 GeV scale. The bulk of the 
Higgs contributions is due to the lightest CP even neutral Hq. As discussed in previous 
sections the dominant Higgs contributions come from the diagrams of Figures 3a,e and 
3f in which a light Higgs is exchanged. These are actually the only sources of Higgs 
contributions in the SM with Hq replaced by the Standard Model Higgs boson. 

The last sector to be discussed is the neutralinos and charginos which in some cases, 
depending on the given inputs, can accomodate light states. Their contributions in 
that case are not necessarily small and is the principal source of deviations from the SM 
predictions. The contributions of the neutralinos and charginos are sensitive to the input 
value for the soft gaugino mass M1/2 being more important for values M1/2 < 100 GeV. 
For such values of the soft gaugino mass and in the unphysical region, < y/s < 
2Mw they are enhanced, due to the development of an anomalous threshold in this 
region, which results to sizable corrections to the magnetic dipole and electric quadrupole 
moments of the W boson 0. 

In Figure 4 and 5 we plot the contributions of the Higgses and of the neutralino - 
chargino sector to Ak^^z, AQr^^z for the most physically interesting case (300, 300, 80) 
GeV and for values of ^/s ranging from GeV to 1 TeV . The region from to 2Mvi/ 
is unphysical since the external have been taken on their mass shell. At s = the 
quantities A/c^, AQ^ are linearly related to the magnetic moment and electric quadrupole 
moments of the W^-boson. 

The structure shown in the Higgs contributions for < 200 GeV is due to the 
lightest of the Higgses. One observes a fast fall off as we increase the energy to values 
above ~ 200 GeV . In the neutralino and chargino sector, and for the /i > case, a 
sharp peak is observed in the unphysical region, -^i < IM-w due to the appearance 
of the anomalous threshold discussed previously and the contributions of this sector is 
substantially enhanced. However such an enhancement does not occur in the physical 
region since their contributions fall rapidly to zero as we depart from the unphysical 
region to values of energies above the two W production threshold. This behaviour is 
clearly seen in Figure 5. The structure observed at energies around 700 Gel^ comes from 
the graph of Figure (3d) and is due to the fact that for these energies we are close to 
thresholds associated with the heavy neutralino states. 

The total contributions to the TGV's both in the MSSM and SM are shown in 
Figures 6 to 8. We display both dispersive and absortive parts of the form factors under 
consideration. One notices that all form factors tend to zero fairly soon with increasing 
the energy reaching their asymptotic values at energies \fs ~ fewTeV in agreement 
with unitarity constraints. 

Our conclusion is that for energies 2M]v < y/s < 200 GeV the MSSM predictions 
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differ in general from those of the SM but they are of the same order of magnitude. 
Deviations from the SM predictions to be detectable require sensitivities reaching the 
per mille level and hence unlikely to be observed at LEP200. If deviations from the SM 
predictions are observed at these energies will be the signal of new underlying dynamics 
which however will not be of supersymmetric nature. At higher energies SM and MSSM 
predictions fall rapidly to zero, due to unitarity, getting smaller by almost an order of 
magnitude already at energies ^/s .5 TeV. As a result, the task of observing deviations 
from the SM which are due to supersymmetry becomes even harder at these energies 
demanding higher experimental accuracies. 



Acknowledgements 

Work supported by EEC Human Capital and Mobihty Program, CHRX - CT93 - 0319. 



17 



A Appendix 

The functions Di[r), D2{r, R) and Q{r) through which the Higgs contributions are ex- 
pressed are defined as follows: 

D,(r) .U\t t 2«^+(-2-r)t3 + (4 + r)t^ 

^' 2 Jo Jo + r{l - 1) - ^t^a{l - a) + te ^ ' 

D2(r,R) = - / dt da ^ ^ -^^ 55 

^ ^ 2 Jo Jo t^ + (-i-r + R)t + r-^t^an~a)+te ^ ' 

Q{r, R) = 2 [' dt [' da ^^(1 - t)a(l - a) 

Jo Jo t^ + [-i-r + R)t + r-^t^a{l-a) + ie ^ ' 



B Appendix 

The Imaginary part of the integral I defined in the main text (see Eq. |53 ) is given by 



Iml=-n [ dt [ da{P2a'^ + Pia + Po) x 
Jo Jo 



40^ 

6{jj^t^a{l-a) + aa + p) (57) 



The roots pi^2 of the argument of the delta function in the expression above are given 

by 



Am = 5I(1-|^)±/(1-|^)'-4|j] . 
The imaginary part Im I vanishes if 

After a straightforward calculation one arrives at the following result 



/m J = -vr ^ dtiP^X^ + Pi A, + Pq) e(A,)e(l - A^) x 

Q{ {st'-a)'-4pt's) 
st^ - (t)2 - 4 pt^s 
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Table Captions 



Table I: MSSM predictions for AA;^^^, ^Q^^^z-, in units of g'^/lSn'^, for three different 
inputs of Aq, mo. Ml/2. Both fj, > and fi < cases are displayed. The energy is 
v^s = 190 GeV. The SM predictions for Higgs masses 50, 100 and 300 GeV respectively 
are also displayed. 

Table II: MSSM predictions for Ak^^z, AQ'y,z, in units of g'^/ldn'^, for three different 
inputs of Aq, mo, Ml/2. Both n > and /x < cases are displayed. The energy is 
^/s = 500 GeV. The SM predictions for Higgs masses 50, 100 and 300 GeV respectively 
are also displayed. 



Figure Captions 

Figure 1: The kinematics of the WWV vertex. 

Figure 2: Fcrmion and sfcrmion contributions to the WWV vertices [figs. a,b]. Neu- 
tralino (Z) and Chargino {C) contributions to WW'j [fig. c] and WWZ [figs. c,d] 
vertices. 

Figure 3: Higgs graphs contributing to the WW'j and WWZ trilinear vertices, ho 
denotes the lightest neutral Higgs boson. 

Figure 4: MSSM Higgs contributions to the dispersive parts of AA;^,, AQ^ (solid lines) 
and Akz, AQz (dashed lines), in units of q'^/IGti'^, as functions of the energy ^/s. The 
inputs are (Ao, mo, M1/2) = (300, 300, 80) Gel^, tan/5 = 2, = 170 GeV. Both > 
and fjL < cases are displayed. The vertical dotted line indicates the position of the 21^ 
production threshold. 

Figure 5: As in Figure 4 for the Neutralino and Chargino contributions. 

Figure 6: MSSM predictions for Ak^^z AQ^ z- The parameters are as in Figure 4. 

Figure 7: MSSM predictions for the absortive parts of Ak^^z, AQ^ ^ ■ The parameters 
are as in Figure 4. 

Figure 8: SM predictions for the dispersive [figs. a,b] and absortive [figs. c,d] parts of 
Ak^^z, AQ^ z fo'^ ^ Standard Model Higgs mass equal to 100 GeV" and mt = 170 GeV. 
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TABLE I 


Ao, mo, Ml/2 


300, 300, 80 


[300, 300, 300] 


(0,0,300) 




tan /3 = 2 


mt = 170 GeV 




= 190 GeV 


/X > 




At < 


AL, 


-1.938 [-1.742] (-1.762) 




-1.733 [-1.767] (-1.791) 


AQ^ 


0.906 [ 0.526] ( 0.538) 




0.300 [ 0.529] ( 0.540) 


Akz 


-2.282 [-2.155] (-2.154) 




-2.016 [-2.146] (-2.140) 


AQz 


-0.354 [ 0.345] ( 0.334) 




-1.193 [ 0.363] ( 0.357) 


SM 


Ak^ = -2.005, -1.735, -2.118 




Akz = -1.350, -2.437, -1.404 


predictions 


AQ^ = .524, .530, .503 




AQz = .507, .533, .481 



TABLE II 


Aq, mo, Ml/2 


300, 300, 80 


[300, 300, 300] 


(0,0,300) 




tan /5 = 2 


= 170 GeV 




= 500 GeV 


/X > 






Ak^ 


-.262 [-.151] (-.191) 




-.310 [-.207] (-.259) 


AQ^ 


0.150 [0.030] (0.056) 




0.146 [0.041] (0.069) 


Akz 


0.121 [0.204] (0.209) 




0.240 [0.191] (0.198) 


AQz 


0.358 [-.407] (-.427) 




0.256 [-.325] (-.352) 


SM 


Ak^ =-.250, -.168, .046 




Akz = .147, .208, -.036 


predictions 


AQ^ = .054, .057, .064 




AQz = .057, .058, .077 
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Figure 1 




Figure 3 



